Nonadiabatic noncyclic geometric phase and ensemble average spectrum of 
conductance in disordered mesoscopic rings with spin-orbit coupling 
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We generalize Yang's theory from the U(l) gauge field to the non-Abelian U(l) x SU(2) spin gauge 
field. Based on this generalization and taking into account the geometric Pancharatnam phase as 
well as an effective Aharonov-Bohm (AB) phase in nonadiabatic noncyclic transport, we calculate 
the ensemble average Fourier spectrum of the conductance in disordered mesoscopic rings connected 
to two leads. Our theory can explain the experimental results reported by Morpurgo et al. [Phys. 
Rev. Lett. 80, f050 (1998)] satisfactorily. In particular, we clarify that the experimentally observed 
splitting, as well as some structure on the sides of the main peak in the ensemble average Fourier 
spectrum, stem from the nonadiabatic noncyclic Pancharatnam phase and the effective AB phase, 
both being dependent on spin-orbit coupling. 

PACS numbers: 73.23.-b, 03.65.Bz, 71.70. Ej 
As is well known, the geometric phase Jl]-^| has manifested itself extensively in physics, particularly in mesoscopic 
systems where quantum interference is extremely important Recently, Morpurgo et al [[Io| reported a novel 

splitting of the main peak (corresponding to the hc/e Aharonov-Bohm (AB) oscillations) in the ensemble average 
Fourier spectrum of the conductance in open mesoscopic rings [ [TT| . The authors conjectured that the observed splitting 
is due to the spin-orbit(SO)-induced Berry's phase |^-^|. It is probably strong experimental evidence showing an 
important effect of the SO geometric phase on quantum transport. Although it was concluded that, in a mesoscopic 
ring possessing the time-reversal symmetry in the absence of AB-flux, the SO-dependent transport can be treated 
formally in the absence of SO-coupling but with an effective magnetic flux ||, it is unclear yet how to calculate the 
value of this flux as well as the ensemble average spectrum of conductance; besides, it is not clear either whether the 
' above conclusion is still valid in the presence of an arbitrary local magnetic field on the ring(i.e., the aforementioned 
time-reversal symmetry is broken), which appears to be the experimental case in fl(i|| . Mal'shukov et al attempted 
to account for the observed splitting of the main peak, but were not quite successful Q. Three aspects of the 
O \ experiment require theoretical explanation. First, the magnitude of the observed splitting is surprisingly large when 
compared with an estimation based on the adiabatic approximation in a clean mesoscopic ring. For a clean ring 
| with radius r subject to a crown-shaped effective magnetic field T5 e ff = {Bocosip r , Bosimp r , B z ) in the cyclindrical 
J> . coordinates, Stern |jj predicted that the inverse B z period becomes (Ai? z ) _1 = 7rr 2 /0 o ± 1/(2_B ) at B » B z , 
0^ ' and the splitting would be ~ 1.2 x 10~ 3 mT~ 1 for an estimated experimental parameter Bq ~ 0.8 T, which is at 
least one order of magnitude less than the observed value. Secondly, the origin of the side structure on the main 
peak needs to be clarified. Lastly, most existing theoretical estimations are crucially based on the adiabatic or cyclic 
condition; however, neither adiabatic nor cyclic evolution is well satisfied in the experiment fipj,[l3). In view of these 
facts, we believe that the nonadiabatic noncyclic geometric phase, essentially similar to the effective flux addressed 
formally in Ref . J5| , could play a crucial role in the system, this is the key point of the present work, being essentially 
different from some existing theoretical analyses (6| J|. However, it is still highly nontrival to evaluate the SO-induced 
geometric phase in the nonadiabatic noncyclic transport and its effect on the ensemble average spectrum. 

It is worth pointing out that the relevant geometric phase detected in the experiment is likely induced by a 
£/(l) x SU (2) sp i n field. In exploring the global geometrical connotations of gauge fields of either the Abelian U (1) type 
or the non-Abelian monopole type, Yang et al showed that the nonintegrable gauge phase factor in the wavefunction 
gives an intrinsic and complete description of the relevant field p3|, hereafter refered to as Yang's theory. In this 
Letter, we first generalize Yang's theory to the non-Abelian £7(1) x SU(2) sp i n electromagnetic field. Using a simple 
one-dimensional(lD) continuum model for a quasi-lD mesoscopic ring, we then analyze carefully the nonintegrable 
phase induced by this field and evaluate its effect on the splitting of the main peak in the ensemble average spectrum 
of the conductance. Remarkably, we find that the splitting as well as the side structure of the main peak observed by 
' Morpurgo et al stem from the nonadiabatic noncyclic geometric phase. 

We consider an electron subject to an electromagnetic field. The corresponding Hamiltonian with U(l) em x 
SU(2) spin gauge symmetry is given by |]| 

H = -^(p+-A-^a) 2 -eA + A1 a° + y(r), (1) 
2m c c 

where /i = g^B/2 with g the gyromagnetic ratio and hb — eh/(2mc) the Bohr magneton. Here A v = (A , A) 
represents a U(l) em electromagnetic potential, and a" = (a , a) = (— a ■ B, er x E/2) is an SU(2) sp i n potential with 
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a denoting the Pauli matrix. V(r) is an arbitrary spin-independent local potential at the point r. The Schrodinger 
equation for the normalized two-component wave function ^(x v ) reads 

iftj^tf (a") = ZT* (sb"). (2) 

By introducing a new wave function p5[ , ^q(x u ) — {/^(x 1 '), where 

U = exp(i— / A l/ dx u )Pexp(i^- / a u dx u ) (3) 
Tic J r Tic J r 

with P the path ordering operator and T an integration curve from a fixed Xq to x", we find that Eq.(^) reduces 
exactly to 

i?i^ (x u ) = H 9 (x v ) (4) 

with 

H„ = m - a §- )6 -i = t^ + V(r) . (5) 

Clearly, U is a continuous local gauge transformation. Under this gauge transformation, the Hamiltonian ([!]) is 
transformed to a Hamiltonian devoid of electromagnetic fields, but with a phase shift in the wave function as seen in 
Eq.(|). In this sense, the gauge factor in Eq. (||) is just the nonintegrable phase in Yang's theory, which can describe 
completely the U(l) x SU(2) sp i n electromagnetic field. For a mesoscopic ring where the phase memory is retained by 
electrons, we may conclude that physical properties of the system in the presence of an electromagnetic field can be 
expressed in terms of the same quantity in the absence of the electromagnetic field, but with a nonintegrable phase 
being taken into account. An important application is related to the SO coupling: any spin-independent transport 
quantity can be expressed in terms of the same quantity in the absence of SO scattering but with an effective magnetic 
flux, a fact which was shown directly by using the transfer matrix method in a tight-binding form for a mesoscopic 
ring possessing the time-reversal symmetry in the absence of AB-flux . In fact, with the help of this generalized 
theory, we are able to study a disordered mesoscopic system subject to an electromagnetic field in a simpler way. 

We now focus on the phase factor first. To capture essential physics of geometric phase in the present quasi-lD 
system, we employ a simple ID model. For a closed path parameterized by arc length s, the total phase factor in 
Eq.(||) is j t = jab +7) where Jab = 2tt(/)/4>o is the usual AB phase with <j> the magnetic flux and </>o — hc/e, and 7 
is the second phase factor in Eq. (g) , which is determined by a Schrodinger- type equation |p|,|l6| 

Here v is a unit vector along the direction of the velocity v = v~v and ds = vdt. Equation (||) describes the evolution 
of the spin state |£) governed by the operator U. The phase associated with Eq.(^) can be further written as 
7 = 7^ + Jab + 7p Hi with the dynamical phase 

l d = ^ j m\° -\m{*))da, (7) 



the effective AB phase 



fli = -f [ * -(v x E)\Z(s))ds, (8) 



and 7 p is the Pancharatnam phase, to be addressed in detail later. Here we emphasize that Jab ^ s a kind of geometric 
phase, though it seems from Eq.(||) as if it were a dynamical phase related to an 'effective magnetic field' — v x E/2c. 
The reason lies in the fact that the two waves propagating in opposite directions in the ring acquire phases with the 
opposite sign for Jab (simply because it depends on the velocity direction v), but the same sign for 7^ p7| . The 
geometrical feature of Jab seems to be ignored in some earlier analyses 1^,^], which appears to be a minor reason 
for the existing discrepancy between theory and experiment. In fact, j A b is just induced by an SU(2) sp i n vector 
potential a, and it is clear from Eq.(|l|) that a plays a role similar to that of the U(l) em vector potential A in the AB 
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effect. As a resuft, it is expected that an effective AB effect can be induced by this SU(2) sp i n vector potential [ p^[ , 
as was also shown by Choi et al 

For a unit vector n — (711,712,^3) = (sindcosip, sin9sirup,cos9) with n £ a unit sphere S 2 , each n corresponds 
to the spin state |£) — {e~ lip / 2 cos{Q /2), e lip ^ 2 sin(9/2)) T via the relation n = (£| <r where T represents matrix 
transposition. The noncyclic Pancharatnam phase accumulated in an evolution of n is found to be |19f | 

1v = -\t n ^ s > (9) 

where dH is an area element on 5 2 , C is a specific closed curve on S 2 , which is along the actual path of n(s) 
plus the shorter geodesic curve from the final point n(s^) = (sinOfCOSipf, sinOfsinipf, cos9f) to the initial point 
n(0) = (sinBiCosipi, sin0isin(pi,cos8i). This Pancharatnam phase can be derived as [|19| 

1 f tf n 1 n 2 -n 2 ni sin(<p f - ipt) 
7p = -o/ TT dt + arctg r — -, (10) 

2 Jo l + n 3 ctgZl-ctgf + cos((p f - (pi) 



where tf is the final time, and n is determined by the equation 



dn 2u ,„ 1 
*=" ^ B -2^ VXE ) Xn ' 

which represents a spin-^ particle moving in an effective magnetic field (B — v x E/2c). This phase is not equal to the 
cyclic Aharonov-Anandan (AA) phase in general [9], but recovers the AA phase jaa — — -| JJ" dt(n\ri2 — n2n\) / {\ + n 3 ) 
for any cyclic evolution with the period r |19| . It is remarkable that the nonintegrable phase in Eq. (|j|) can be evaluated 
by simply computing Eqs.(|7|), (||), and ([l0); while it is hard to calculate the value of the effective flux addressed formally 
in Ref. ||, particularly in the presence of an arbitrary local magnetic field. 

At this stage, as in the experiment fl(i[| , we study a disordered ring with the Rashba SO-interaction(equivalent to 
an internal electric field E = Ee z ), subject to a local magnetic field B = B z e z and a magnetic flux <fi — nr 2 B z . The 
Hamiltonian, which is in the form of Eq.(@), becomes pjl3|| 

H = Tiuj r {-i- h -(a x cosip r + a y sin<p r )] 2 - [iB z a z + V(ip r ), (12) 

Oip r 00 2 

where u> r = h/ (2mr 2 ), ip r is the polar angle, and the normalized electric field strength r\ = [isEr/ch = Imnr with 
the SO coefficient Ti 2 k. 

To account for the experimental results naturally, we investigate the electronic transmission across a disordered ring 
connected to external current leads, schematically illustrated in Fig. 3 in Ref. |20| . In such a system, the electronic 
transmission is significantly affected by the nonintegrable phase. Using the method originally proposed by Biittiker 



et al j 20 j and our generalization of Yang's theory, the transmission coefficient across the ring is found to be 

(h - n -nT. r. , , , , , , 

b 2 \ -a 1 J V -a 1 I + 1 \ -1 



T — — 

9 ~ 6 4 



(13) 



where T+ and T_ are the transfer matrices of the upper and lower branches of the ring, 1 is the unit matrix, 
a = ±(Vl - 2e - l)/2, and b = ±(vT - 2e + l)/2 with < e < 1/2. A7 = 7^3 + 7^/(n(0)) + 7 P (n(0)) represents 
the nonadiabatic noncyclic geometrical phase accumulated in the evolution when electron(with the initial spin-state 
n(0)) moves one cycle in the clockwise sense pi] ]. For a beam of electron waves with Fermi wave vector kf, the rate 
for electrons to traverse one round in the ring is w/ = hkf /(mr) for ballistic motion, but is estimated approximately 
to be ujd — lujf/{2'Kr) for weak diffusive motion [p2| , where I is the electron mean free path. This rate can be 
regarded as the angular frequency of the otherwise rotating magnetic field felt by the electron spin p2| , which is 
given by B e yj(t) = (B^cosLOf^t, BQ^sinujf^t, B z ) with B^ d — —rjfiuj f^/2/i. Then from the equation dn(t)/dt = 
— (2{j,/fr)Ti e f f(t) x n(t), n(t) is derived exactly as 

(cosujf^dt —sinojf.dt 
sinujf.dt cosujf^dt 
1 

sin 2 x + cos 2 xcosu> s t cosxsinui s t ^sin2\(l — cosuj s t) 

I —cosxsinu s t cosuo s t sinxsinuj s t | n T (0), (14) 

^sm2%(l — cosu> s t) —sinxsincu s t cos 2 x + sin 2 xcosu s t 
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where uj s = y/tOQ + {^f.d + ^i) 2 and \ = arc ^5[ w o/(^/,d + w i)] with loq — 2/iB^ d /U and wi = 2fj,B z /T%. On the other 
hand, we can rewrite Eq.(^) clearly as 

,, nojr,! r 27r / u f. d 

7 7/(n(0)) = — '-M / sin9cos{u f , d t - <p)dt (15) 



with 9 = arctg{^J n 2 + n 2 /n 3 ) and (/? = arctg(n 2 /ni). Substituting Eq.(|l4|) into Eqs.([To|) and (|l5|), the nonadiabatic 
noncyclic phases and 7 P can be computed, at least numerically. 

For simplicity, but without loss of generality, we compute T + and T_ in a generalized Kronig-Penny ring consisting 
of N = N + + N- uniformly spaced 8— function barriers with random strengths, where N + (N-) is the number of 
barriers on the upper (lower) branch. The Hamiltonian for the system in the absence of electromagnetic fields reads 
Hq = — (h 2 /2m)d 2 /dx 2 + ^(- T ? — n^ao), where A„ s is a potential strength parameter, ao is the lattice 

spacing, and <; = +(— ) represents the upper (lower) branch. The spinless electron wave function in the regions where 
no potentials are present may be written as ip^x^) = A n( e lk f x,; -I- B n<; e~ f x,! . The coefficients A n and B„ s across 

site x<; = n c Oo are related through the matrix M„ s , 

£c»+u« ) = m % ( i n 

-D(n+l) s - / V ti n 



with 



/ 1 _ ^ »Vn t -2ik ( n c a n 

M - 2k f 2k f 

1U n, ~ iV n( 2ik f n c a n i , 

\ 2k f C T 2fc/ 



where V^ s = 2m\ n ^/h 2 is assumed to be distributed uniformly in an interval [— w/2, w/2\. Then, one can find 

~ / p ik f N,a n \ — 



n^ — 1 



Note that T ? can be further simplified to a 2 x 2 matrix |23(|. Substituting the simplified Eq.(16) and A7 into Eq.(|13| 



we are able to calculate the transmission coefficient T g 



For comparison with the experimental observation [ jro| ], we plot in Fig.l the calculated ensemble average Fourier 
spectrum of the conductance for unpolarized electrons , which is defined as 

(\G(u)\) = (\ I**" e™ B *G{B z )dB z \), (17) 



where G{B Z ) — (e 2 /h) Yl±a(o) T g (B z ) with T g the average on the initial spin-orientation (24|, ( ) represents the 
ensemble average. Reasonable parameters in the calculation are determined as follows. As in the experiment, B m = 
0.35 T, v f ~ 3.0 x 10 5 m/s, g ~ 14, r ~ 1.05jum (with A ? = 4200 and a ~ § X 10~ 9 m), which leads to the 
period in a magnetic field ~ 1.2 mT; w ~ 0.267&/, which corresponds to the mean free path I = 96k 2 ao/w 2 ~ 
l.O^m pq ). The dimensionless coefficient rj ~ 3.5, which corresponds to the experimentally reported SO coefficient 
h 2 K ~ 5.5 x 10~ 10 eVcm |p0| . Finally, it is typical to consider the case e = 0.25. It is worth emphasizing that the 
essential feature of Fig.l is sensitive mainly to the SO coupling parameter rj: no clear splitting is present in the main 
peak if r\ is smaller than about 1.5. This implies that the SO-interaction plays a crucial role in the splitting. To 
clarify the origin of the structure of the main peak, we plot it under both adiabatic and nonadiabatic conditions. 
The former case is shown in the inset of Fig.l, where 7 P = ircos9i(l — cos\ a ) and = —r]Trcos9iSinXa with 

Xa = arctg(B^/B z ). As (B z /Bq) is no longer small, it is very difficult to have a clear analytical understanding of the 
influences of 7 P and 7^3 on the splitting of the peak. From the inset of Fig.l, we can see that under the adiabatic 
approximation a somewhat splitting of the main Fourier peak is present only if we include the effective AB phase; 
however, this splitting feature is obviously not in good agreement with the experimental result. After careful analysis, 
we understand that the effect of the adiabatic 7 p -phase is too weak to play an important role in causing clearly 
observable splitting. More remarkably, if ever we take into account both the Pancharatnam phase and the effective 
AB phase in the nonadiabatic noncyclic case, as shown by the solid line in the main panel of Fig.l, our theoretical 
result is in excellent agreement with the experimental observation, especially for the splitting and the side structure 
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(two small peaks) of the main peak (see Fig. 5 in Ref. |lCj]). From Fig.l, we can also see that the Pancharatnam phase 
plays a key role in the main splitting, while the two small side-peaks are closely related to the effective AB phase, 
ft is therefore clarified for the first time that the experimentally observed splitting of the main peak in the ensemble 
average Fourier spectrum stems from the nonadiabatic noncyclic Pancharatnam phase and the effective AB phase, 
both being dependent on the SO coupling. 

Finally, we remark that the radius of the ring ~ the mean free path (~ lfim) and the transport is in the weakly 
diffusive regime in the experiment. In the present quasi-lD ring, a multi-channel effect, albert weak and secondary, 
may exhibit in the ensemble average spectrum of the conductance (e.g., broadening and smearing of the peak splitting), 
which has been ignored in the present work and may deserve for further detailed study in future. Nevertheless, the 
effect would not affect the present conclusion regarding the splitting qualitatively because the non-adiabaticity of 
geometric phase is unlikely changed significantly in the present weak diffusive ring. 

* To whom correspondence should be addressed. E-mail: zwang@hkucc.hku.hk 
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Figure Caption 

Fig.l The peak of the ensemble average Fourier spectrum of the conductance in nonadiabatic noncyclic cases. The 
inset shows the corresponding curves under the adiabatic condition. 
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